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Three-dimensional topological photonic crystal with a single surface Dirac cone 
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A single Dirac cone on the surface is the hallmark of three-dimensional (3D) topological insulators, where 
the double degeneracy at the Dirac point is protected by time-reversal symmetry and the spin-splitting away 
from the point is provided by the spin-orbital coupling. Here we predict a single Dirac-cone surface state in a 
3D photonic crystal, where the degeneracy at the Dirac point is protected by a nonsymmorphic glide reflection 
and the linear splitting away from it is enabled by breaking time-reversal symmetry. Such a gapless surface 
state is fully robust against random disorder of any type. This bosonic topological band structure is achieved by 
applying alternating magnetization to gap out the 3D “generalized Dirac points” discovered in the bulk of our 
crystal. The Z 2 bulk invariant is characterized through the evolution of Wannier centers. Our proposal—readily 
realizable using ferrimagnetic materials at microwave frequencies—can also be regarded as the photonic analog 
of topological crystalline insulators, providing the first 3D bosonic symmetry-protected topological system. 


Topology of electron wavefunctions was first introduced to 
characterize the quantum Hall states in two dimensions (2D) 
discovered in 1980 m. Over the past decade, it has been 
recognized that symmetry plays a crucial role in the classi¬ 
fication of topological phases, leading to the broad notion 
of symmetry-protected topological phases As a pri¬ 

mary example, topological insulators lEHII are distinguished 
from normal insulators in the presence of time-reversal sym¬ 
metry (T). A 3D topological insulator exhibits an odd num¬ 
ber of protected surface Dirac cones, a unique property that 
cannot be realized in any 2D systems. Importantly, the exis¬ 
tence of topological insulators requires Kramers’ degeneracy 
in spin-orbit coupled electronic materials; this forbids any di¬ 
rect analogue in boson systems IF]. In this report, we discover 
a 3D topological photonic crystal phase hosting a single sur¬ 
face Dirac cone, which is protected by a crystal symmetry lEl- 
iia — the nonsymmorphic glide reflection rather than T. Our 
finding expands the scope of 3D topological materials from 
fermions to bosons. 

Unlike in Fermi systems, achieving a single Dirac cone in 
boson systems requires T breaking. This is because the T op¬ 
erator acts differently on bosons and fermions; Tf = — I for 
fermions with half-integer spins and i for bosons with 

integer spins. As a result, Tb is not compatible with the Hamil¬ 
tonian of a single Dirac cone, while 7/ is (see Supplemen¬ 
tary Information). Instead of Tf, the Dirac point-degeneracy 
in our photonic crystal is protected by a glide reflection ini- 
na , which ensures a band-crossing on high-symmetry lines. 
This crystal-symmetry-protected topological photonic crystal 
can be regarded as a bosonic analog of the recently discovered 
topological crystalline insulators in electronic systems 13 [Bj- 

M- 

Our starting point is a photonic crystal having a body- 
centered-cubic (bcc) unit cell which contains four identical di¬ 
electric rods, illustrated with different colors for clarity in Fig. 

. This crystal belongs to the nonsymmorphic space group of 
No. 230 (Ia3d) that contains glide reflections and inversion. 
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Interestingly, such a triply-periodic structure self-assembles 
as disclination-line networks in the first blue phase of liquid 
crystals denoted as BPI. Here, the dielectric constant (e) 
of the rods is 11 and radius is 0.13a where a is the length of 
the cubic cell. In Fig. [T]\ , the photonic band structure of BPI 
shows a four-fold degenerate point at the P momentum, dis¬ 
persing linearly in all directions of the 3D momentum space. 
Unlike a 3D Dirac point ll20l 1211 — a four-fold degeneracy 
point which splits into two sets of doublet bands along any di¬ 
rection — our four-fold degeneracy here splits into four bands 
along a generic direction. In Fig. [^, this splitting is not ob¬ 
vious, since most dispersions still remain doubly-degenerate 
along high-symmetry momentum lines. However, it is clear 
that the 3rd and 4th bands split along P — F and the 1st and 
2nd bands split along P — H. We name this type of degeneracy 
a 3D “generalized Dirac point” (GDP). We note that there are 
two non-equivalent P points in the 3D bcc Brillouin zone (BZ) 
related by inversion. Interestingly, the two GDPs (at ±P) are 
the only band-touchings between band 1,2 and 3,4. When the 
space group is perturbed, the GDPs could turn into line nodes, 
Weyl points fT2\ 1^ or open bandgaps. Detailed studies of 
GDPs will be presented in another paper. 

In symmorphic space groups, where the point groups de¬ 
couple from lattice translations, the highest dimension of 
group representation is three. The four-fold band degen¬ 
eracies of the GDPs are hence the consequence of the non¬ 
symmorphic symmetries of glide reflections and screw rota¬ 
tions in BPI. A non-symmorphic symmetry is in general com¬ 
posed of a point group (mirror or rotation) followed by a frac¬ 
tional lattice translation, where neither of the two is a symme¬ 
try of the system. The important feature of a non-symmorphic 
space group is the extra band degeneracies at the BZ bound¬ 
aries EMI. Since the screw rotations cannot be preserved 
on a planar surface, we focus on the glide reflections in or¬ 
der to obtain protected surface states. Shown in Fig. [Tp, the 
(001) surface has two glide reflections Gx = {Mrijx— jy} 
and Gy — {My\^x}. The top view schematic illustrates the re¬ 
lations between the four rods under the two glide reflections. 
The (001) surface BZ is plotted in Fig. [Tp. 

A glide reflection ensures a linear point-degeneracy along 
each glide-reflection-invariant momentum line. To see this, 
we study the Bloch states on the Gy-invariant lines ofX'—X 
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FIG. 1. Bulk band structures of the BPI photonic crystal. A) The GDP is between the first four bands. B) The cubic unit cell of length a 
consisting of four identical dielectric rods oriented along the bcc lattice vectors of (111) (red), (Ill) (yellow), (ill) (blue) and (111) (green). 
The rods go through (0,0,0)a, (0,0.5,0)a. (0.5,0,0)a and (0,0,0.5)a respectively. There are two glide reflection planes (Gx and Gy) in the 
structure, invariant on the (001) surface. The top-view schematic illustrates the relations between the rods under operations of Gx and Gy. C) 
The bcc BZ and its projection onto the (001) surface BZ. The transparent red box outlines the volume in the bulk BZ that projects to half of 
the surface BZ. D) Bulk band structure showing that the GDP opens when magnetization is applied on the rods without breaking Gy. 


and M’ —M dashed in the (001) surface BZ on the right of 
Fig. 1^. A Bloch state with momentum {kx,ky) is mapped 
to another state with momentum {kx,—ky) under Gy, so for 
any state along these two lines with ky = 0 and ky = Tija, its 
momentum is invariant under Gy. This means the eigenval¬ 
ues of Gy igy) are good quantum numbers for the Bloch states 
on these two lines. Since Gy'^ = {l|ax}, gy{kx) = 

[gy{kx) = which is A:x-dependent. The two branches of 

glide-reflection eigenvalues always differ by a minus sign and 
they evolve into each other after a 27r transportation along the 
Gy-invariant lines due to the fact that gy{k) = —gy{k+ ^). As 
a result, the corresponding wavefunctions of the two branches 
have the same winding as their eigenvalues—a unique prop¬ 
erty of the half-lattice translation in glide reflections. Con¬ 
sequently, the two frequency eigenvalues of the two Bloch 
modes also switch values after transporting a period along the 
invariant momentum lines, illustrated in Fig. [^. Assume the 
frequencies of the two modes are co^ and co at an arbitrary 
kx point (say kxU = 0). The frequency dispersions switch their 
values at ka = 2n. This switch ensures a crossing point (red 
dot) onX' — X and M' —M respectively. We argue that these 
two protected double-degeneracies give a Z 2 classification of 
the surface states ifT^ . Illustrated in the middle of|^, there 
are two topologically in-equivalent ways for these two point- 
degeneracies to connect. The gapless connection is a signature 
of the topologically nontrivial surface states protected by Gy. 

We now break T in the BPI photonic crystal to open the 
bulk bandgap without breaking the Gy. Shown in Fig. [IP’ 
the GDP at P point lift up into a bandgap when we apply al¬ 
ternating magnetization on the rods along z. These magne¬ 
tization induces off-diagonal imaginary parts in the dielectric 
constant (e) of materials with gyroelectric response j^ . (Fer- 
rimagnetic materials of gyromagnetic response 12911 give the 
same results in Supplementary Information). Here jJ. = 1 and 


/ £// (f 0 \ 

e = ~K Zjj 0 , where = 11, — \k\^ = e'^ ||22l 

VO 0 £zzj 

and Jf is a non-zero imaginary number when the magnetiza¬ 
tion (M^) is present. In Fig. EP’ K = —10/,—5/,+5/,+10/ for 
the red, yellow, blue and green rods respectively. This config¬ 
uration preserves Gy, because magnetization (magnetic held) 
flips sign under a mirror (reflection) operation. The 2D plane 
group of the resulting (001) surface is pg. 

The (001) surface state, plotted in Fig. |2fi, has a single 
Dirac cone at point L on the M' —M line, consistent with the 
glide-reflection degeneracy in Fig. |g\. By varying the mag¬ 
netization or rod radius without breaking Gy, the Dirac point 
L moves along the Gy invariant line M' —M. This single Dirac 
cone at L is connected gaplessly with the bulk bands across 
the band gap. In Fig. W: we restore Gx to coexist with Gy by 
doubling the magnetization amplitude of the green and yel¬ 
low rods (I Jf I from 5 to 10). The surface plane group becomes 
p2gg. Due to this extra glide-reflection plane through Y point, 
the surface Dirac cone is then pinned at Y on M' — M. If 
we break both glide symmetries by de-magnetizing the yel¬ 
low rod, both glide planes of Gx and Gy are broken and the 
surface plane group reduces to pi. The surface Dirac cone is 
now gapped as shown in Fig. HP- This demonstrates that the 
gapless surface states are indeed protected by the glide reflec¬ 
tion. 

The principle of bulk-edge correspondence says that the 
surface state is a holographic representation of the bulk topol¬ 
ogy. We demonstrate this correspondence between the sur¬ 
face states in Fig. |^and the “hybrid Wannier centers” OOl of 
the bulk bands below the bandgap computed in Fi g. This 
approach is also known as the Wilson loops EUlla. The 
hybrid Wannier function of each band is a spatially-localized 
wavefunction along z, obtained from Fourier-transforming the 
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FIG. 2. The (001) surface states after breaking T. The surface 
and bulk dispersions are plotted in red and blue colors, respectively. 

A) Illustration of the two crossing points (double-degeneracies) on 
the Gv-invariant lines of M' — M and X' —X (dotted) in the surface 
BZ. Consequently, the surface states can have a gapless connectivity 
which is topologically nontrivial. M' and M are the same point in the 
BZ, so are X' and X. B) Single Dirac cone at L movable along the 
M' — M line protected by G,.. C) Single Dirac cone in B) pinned at 
Y point due to the co-existence of Gx and Gy. D) Surface states in 

B) and C) gapped by breaking the glide-reflection symmetries. For 
all above surface calculations, a perfect metallic boundary is placed 
on top, at the center of the cubic cell, on the (001) surface with an 
air gap spacing of 0.5a from the photonic crystal surface. The 3D 
plots of the surface dispersions are plotted on the right to show the 
behavior of the surface Dirac cone. The 3D surface plots are centered 
at the Y point, with a span of Q.2n/a in kx and 0.1;r/a in ky and a 
normalized frequency range between 0.41 and 0.43. 


Bloch wavefunctions with respect to while keeping the 
other two surface momenta. The z-position expectation values 
of the hybrid Wannier wavefuntions, i.e., the hybrid Wannier 
centers, are equivalent to the Berry phases of the bulk bands 
below the gap along a loop in z in the bulk BZ. In our bcc lat¬ 
tice, this non-contractable loop (of length Anja) is the vector 
connecting H and —H in Fig. [Tp. This hybrid Wannier center 
is well defined up to a lattice period of a/2 in z, and similarly, 



FIG. 3. Hybrid Wannier centers in the surface BZ indicating bulk 
topologies and the connections of the surface states. (A) The gap¬ 
less hybrid Wannier centers corresponds to the non-trivial surface 
states ipg) in Fig.|^. (B) The gapped hybrid Wannier centers corre¬ 
sponds to the trivial surface states (pi) in Fig. [^. The hybrid Wan¬ 
nier centers corresponds to the p2gg surface in Fig. is plotted in 
the Supplementary Information. 


the Berry phase has a 2% phase ambiguity. The calculations 
of the Berry phases are detailed in the Supplementary Infor¬ 
mation. 

In the surface BZ, a gapless spectrum of Wannier cen¬ 
ters (Berry phases) indicates a non-trivial bulk topology and 
a gapless surface state. In contrast, a gapped spectrum repre¬ 
sents a trivial bulk topology and the absence of gapless surface 
states. This can be understood by the following intuitive ar¬ 
guments. If there is a full gap in the spectrum of Wannier 
centers, then there is a certain position in z where no state is 
localized. Terminating the bulk at that plane results in a sur¬ 
face without surface states — the trivial smTace states. On 
the other hand, if the Wannier center plot is gapless, then for 
any surface termination there must be a localized surface state 
at some surface momentum. The surface is hence gapless for 
terminations at arbitrary z — a nontrivial surface. In Fig. 
we plot the Wannier centers of the two lowest bands along the 
closed loop of X' —X —M — M' —X' in the surface BZ. Fig. 
1^ depicts the hybrid Wannier centers calculated for the bulk 
bands in Fig. [IP^ whose surface state is shown in Fig. IIP- 
Similarly, the hybrid Wannier centers in Fig. correspond 
to the surface states shown in Fig. |2p. The Wannier centers 
are gapless in Fig. 13^. consistent with the existence of the 
gapless single surface Dirac cone in Fig. |^. In comparison, 
the Wannier centers in Fig. is gapped, also consistent with 
the absence of topological surface states in Fig. HP- These 
data confirm the bulk-edge correspondence that the Wannier 
centers for all bulk bands below the bandgap is homotopic to 
the surface band structure of a semi-infinite system with one 
open surface. 

Single-Dirac-cone surface states are fully robust and do not 
localize under arbitrary random disorders on the surface. This 
has been discussed in 3D topological insulators where the sur¬ 
face states remain delocalized under random impurities of any 
type Il3^l34ll . In our case, although individual defects break 
the glide reflection, their ensemble average do not. Intuitively, 
if one local disorder generates a positive Dirac mass term 
within a region on the surface, there must be a neighboring re¬ 
gion where the mass term is negative. A chiral edge mode ex¬ 
ists along the edge between two regions with opposite masses, 
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similar to the photonic one-way edge states Il28ll29l analogous 
to the quantum Hall effect. In the presence of strong disorder, 
these chiral edge modes percolate the surface and the surface 
states remain delocalized. The surface with a strong random 
disorder can be mapped to the electronic states at the critical 
point of a quantum Hall plateau transition, where chiral edge 
modes between regions of different Landau-level filling fac¬ 
tors percolate. The transmission rate of light on the surface 
hence exhibits the universal scaling laws in the universality 
class of the quantum Hall plateau transitions lfT2l iTSll . Free 
from any interaction, this single-Dirac-cone surface state is an 
ideal platform for studying the critical phenomena of “metal- 
insulator” transitions in Dirac systems |[33]|36l. 

In 2D photonic crystals, topological band structrues pro¬ 
tected by e — /z symmetry llJTl f38l have been studied. How¬ 
ever, symmetries in constitutive relations are difficult to main¬ 
tain over a wide frequency bandwidth. Another 2D example 
discusses the bulk topology of Q rotation ||3^ . Unfortunately, 
six-fold rotation cannot be preserved on the ID edge and can¬ 
not protect edge states. In contrast, our glide reflection can be 
maintained for all materials at all frequencies with protected 
surface states. 

Experimentally, the T-breaking BPI photonic crystals can 
be readily realized at microwave frequencies by assembling 


ferrimagnetic rods ll29l l40ll with internal remnant magnetiza¬ 
tion without the need for external magnetic fields. These ma¬ 
terials are commercially available, such as yttrium iron garnet. 
Towards optical frequencies, T-breaking could potentially be 
implemented through dynamic Floquet modulations BTl l42l . 
In addition, our approach for photons can directly be used for 
phonons where T-breaking can be achieved by spinning the 
rods ll43]l . 

This work demonstrates that symmetry-protected 3D topo¬ 
logical bandgaps supporting disorder-immune surface states 
can be obtained in bosonic systems. Spatial symmetries ||9] 
[mmiiiiiia other than the glide reflection are to be stud¬ 
ied in the rich context of 230 space groups and 1651 magnetic 
groups for any bosonic particles. 
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SUPPLEMENTARY INEORMATION 


A. Compatibility between time-reversal symmetries and the 
single-Dlrac-cone surface state 

The key ingredient in achieving bosonic single Dirac sur¬ 
face states is the breaking of T, which we prove by contra¬ 
diction. The T operator differs fundamentally for particles 
with different spins ; Tf = — I for fermions with half-integer 
spins while — I for bosons with integer spins. Up to a 
choice of basis, the anti-unitary T operator can always be ex¬ 
pressed asTf = and % = |k^-k for fermions and 

bosons, respectively. Here are the Pauli matrices acting 
on the two-component wavefunctions of a single Dirac cone. 
T flips the sign of momentum (k) and K is the complex con¬ 
jugation. Let us consider the surface (say the xy plane) of 
a 3D system with a bulk gap. In the presence of T, a single 
Dirac cone can only appear at a T-invariant momentum, in the 
vicinity of which the two-band Dirac Hamiltonian is denoted 
by Hsd{^) = +kyG 2 . Here kj,y are the two surface mo¬ 
menta and (7i 2 are two linearly independent Pauli matrices. 
T-invariance implies T and //yo(k) commute, requiring the 
existence of at least two Pauli matrices that anti-commute with 
T. For fermions, this is satisfied since all three Pauli matrices 
anti-commute with 7/. The anti-commutation relations forbid 
any mass term ( 03 ) in //s£)(k), justifying the T-protected sin¬ 
gle Dirac surface states found in topological insulators. For 
bosons, however, Oj. is the only Pauli matrix anti-commuting 
with %■ So % is not compatible with Hsd- Hence T has to be 
broken to linearly split the two bands in all surface directions 
away from the Dirac point, in order to form a single surface 
Dirac cone in photonic crystals. 


B. Eerrimagnetic materials 


Here we break T with the gyromagnetic material instead 
of the gyroelectric material in the main text. For example, 
yttrium iron garnet (YIG) has strong gyromagnetic responses 
at microwave frequencies. The permittivity and permeabil- 

/ M// V 0 

ity of the YIG crystal can bee = ll, /r= —v /iy 0 

V 0 0 )io 

where v is a non-zero imaginary number when the magneti¬ 
zation (M^) is present and jio = 1. In Fig. SI = 1.5 and 


V = —1.21,— 1.21,-1-1.2/,-1-1.21 for the red, yellow, blue and 
green rods respectively. The bulk bandgap opens and the sur¬ 
face state has the same nontrivial connectivity as that is shown 
in Fig. 15 :. These calculations were performed using a modi¬ 
fied version of the MIT Photonic Bands 


C. Discussions of hybrid Wannier-center spectra 

Using the notation in Ref. lf30l . the hybrid Wannier centers 
Zn{k±) and the equivalent Berry phases (/)„(kx) is related by 
Zn(k±) = ^(j)„ik±). Here n labels the band, k± means the 
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“Pinned” single surface Dirac cone (p2gg) 



X Y X’ M’ Y M X 


FIG. SI. A) Bulk band structure showing the gap openning of the 
GDP. B) Gapless Z 2 surface states of single Dirac cones atX=X’ and 
Y=Y'. A perfect metallic boundary is placed from top, at the center 
of the cubic cell, on the (001) surface with an air gap spacing of 0.4fl 
from the photonic crystal surface. 
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FIG. S2. The gapless hybrid Wannier centers corresponds to the non¬ 
trivial surface states (p2gg) in Fig. 


wavevector in the surface BZ perpendicular to k^- c is the 
real-space period in the direction of the surface, which is a /2 
in our bcc lattice for ( 001 ) surface. 

Hybrid Wannier centers [zn{k±)] are properties of the bulk. 
In addition to the symmetries on the surface, z„ {k± ) also have 
z to —z symmetries in the bulk. In our system, this z-symmetry 
is inversion, so Znik±) = —Zn{—k_i) in all the plots in Fig. 


and S2 At inversion-invariant momenta {X,M,r, Y) where 
k = —k up to a reciprocal vector, £z„(kx) = 0 or ;r. If the 


Wannier centers have a single crossing in the closed loop of 
X' — T — X or M — Y — M', the crossing point must locate at 
one of the inversion-invariant points. The above arguments 
are all consistent with the three plots of Fig. [^,B and Fig. 

m 


Fig. 1^ plots the hybrid Wannier centers in the surface BZ 
of surface symmetry p2gg. Due to the high symmetry, the two 
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Wannier centers are completely degenerate on the X' — Y — X 
line at the phase value of n. The 2D vertical plane in the bulk 
BZ containing X' — Y — X is the only plane on which every k 
point is invariant under Gy. Consequently, the lowest two bulk 
bands on this vertical plane can be separately labeled by the 
two Gy eigenvalues of g^. On the other hand, the multiplica¬ 
tion of Gx and inversion (P) also transforms the Bloch states, 
within this vertical plane, from {kx,Q,k^) into {kx,Q,—k^). In 
addition, [GxP, Gy] = 0. This commutation relation means that 
these two operators share the same eigenstates on the plane, 
so that we can transforming the states within each separate 
branch of the two bulk bands labeled by g^. Since GxP trans¬ 
forms k^ into —k^ in the plane. This z-symmetry requires 
z(kx) = -z(kx),i.e. Berry phases of 0 or ;r for both branches 


of the two bulk bands, for kx on the X' — Y — X line. 

D. Calculation of Berry phases (hybrid Wannier centers) 

The multi-band non-Abelian Berry phases are calculated 
through the linking matrices of the Bloch wavefunc- 

k k I 

tions Umk between the neighbouring two k points. Mmn = 
(Mmk|Mn(k+Ak))- We multiply the linking matrices to be the 

Wilson loop W(kx) = along the closed loop 

in the BZ- a parallel transport cycle. The Wilson loop 
eigenvalues are A„(kx) and the Berry phases are 0„(kx) = 
Im[logZ„(kx)]. 

The key of this calculation is fixing the periodic gauge at 
the two end points (ko and kiast) differ by a reciprocal vector 
G in the bulk BZ. We set M(kiast) = M(ko -fG) = e^'*' ''M(ko). 
For other k points in the loop, gauge fixing is not required. 



